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The effect of spatiai modulation of the single-site anisotropy D on the ground state of the 
5 = 1 Heisenberg chains is investigated. In the case of period 2 modulation, it is found that the 
phase diagram contains the Haldane phase, large- _D phase, Neel phase of udud-type and uOdO- 
type. It is shown that the hidden antiferromagnetic order in the Haldane phase compatible 
with the spatial modulation of D-term get frozen resulting in the emergence of various types 
of Neel orders. The investigation of the model with longer period D-modulation also confirms 
this picture. 
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1. Introduction 

Antiferromagnetic quantum spin chains have been con- 
tinually the subject of experimental and theoretical re- 
search in the past two decades. Among them, the Hal- 
dane ground state of the 5 = 1 antiferromagnetic Heisen- 
berg chain with a gap to the first excited state 1 has 
been extensively studied by many authors. This state is 
characterized by the hidden antiferromagnetic string or- 
der accompanied by the Zi x Zi symmetry breakdown. 
The easy plane single-site anisotropy D{> 0) destroys 
the Haldane ground state leading to the large-D state 
without specific order. On the other hand, the easy axis 
single-site anisotropy (D < 0) drives the Haldane state 
into the Neel state. 2-4 

In this context, it is an interesting issue to investigate 
how the Haldane phase is modified if the easy-axis and 
easy-plane Z?-terms coexist in a single chain. As a sim- 
plest example of such competition, we first investigate 
the ground state phase diagram of the 5 = 1 antiferro- 
magnetic Heisenberg chain with alternating D-term, 
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where D+ = Do+SD, D- = Dq — 5D and Si is the spin-1 
operator on the i-th. site. The parameters Do and 5D rep- 
resent uniform and alternating components of uniaxial 
single-ion anisotropy, respectively. The periodic bound- 
ary condition is assumed unless specifically mentioned. 

The most remarkable conclusion on the ground state of 
this model obtained in the present work is that the uOdO- 
type long period Neel state emerges for large SD although 
this type of order is not compatible with the short range 
order in the Haldane phase at D + = D- = 0. Here u, 
and d stand for the single-site states \Sf) with Sf = 1,0 



and — 1, respectively. This implies that the single-site 
anisotropy can reverse the sign of spin-spin correlation in 
the corresponding isotropic model. This phenomenon can 
be understood as a selection of the specific antiferromag- 
netic order compatible with the D- alternation among the 
components of the Haldane state which is the superposi- 
tion of various states with the hidden antiferromagnetic 
order. We also confirm this idea by investigating the sim- 
ilar model with long period D-modulation. 

This paper is organized as follows. In the next section, 
the effective Hamiltonians in various limiting cases are 
derived and the overall feature of the phase diagram is 
explained. The numerical exact diagonalization results 
and the quantitative phase diagram are presented in §3 
. The model with long period D-modulation is discussed 
in §4. In the final section, we summarize our results. 

2. Effective Hamiltonians 

Let us first consider the various limiting cases of the 
Hamiltonian (1) in which the ground states can be de- 
termined by the effective Hamiltonians derived perturba- 
tionally. In the simplest limiting case of D+,D- >> J, 
all spins are confined to the state |0) and the ground 
state is obviously the large-D phase. 

For D + , D- << — J, all spins are allowed to be in the 
states |±1). The effective coupling between these spins 
are given by 
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up to the lowest order in J. Therefore we expect the 
udud-type Neel phase. 

In the limit of D + >> J, Z?_ << — J, the spins on the 
odd-th sites are confined into the state |0) and those on 
the even-th sites can take |±1) states. The latter spins 
are coupled via the Ising type antiferromagnetic coupling 
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within the third order perturbation in J as, 
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although the spins S^i and 62(1+1) axe ferromagnetically 
correlated in the absence of D-terms. The origin of the 
antifcrromagnetic interaction is as follows. Let us denote 
the spin state of 3 successive sites of the original Hamil- 
tonian (1) by 5|» •S'fi+i £2(1+1))' ^ n tnc absence of the 
exchange term J, the ground states are ±1 ±1} which 
are 4- fold degenerate. Due to the spin flip term in (1), the 
1 1 1 ) state is mixed with the 1 1 1 0) and 1 1} states 
while the |1 —1} state is mixed with the |1 —1 0) and 
|0 1 —1) states. Therefore the energy of the 1 1 1 ) state 
is raised by the nearest neighbour pair of up spins in 
the intermediate state while the energy of the |1 —1) 
state is lowered by the nearest neighbour pair of up and 
down spins in the intermediate state. This means that 
the effective interaction between Sfj and ^(i+i) 1S anti- 
ferromagnetic. It should be noted that the spin nip terms 
in the effective Hamiltonian appear only in the 4-th or- 
der perturbation in J. Therefore we expect the uOdO-type 
Neel phase in this regime. 

This can be also interpreted using the concept of the 
hidden antiferromagnetic string order in the following 
way. In the absence of the D-terms, the ground state has 
a hidden string order which implies that the spins with 
|±1) are arranged antifcrromagnetically if the sites with 
|0) are skipped. 2,3 The position of the sites with |±1) 
and |0) strongly fluctuate quantum mechanically and this 
antifcrromagnetic order remains hidden because it is im- 
possible to observe the correlation between only the sites 
with |±1) experimentally. In the presence of strong D- 
terms, only the states consistent with the constraint set 
by the D-terms survive among all states with hidden or- 
der. For D + , D- « — J, no |0) sites are allowed, so that 
the hidden order is frozen into the explicit udud order. 
For D + » J and D- << J, the odd-th site must be in 
the state |0) and the even-th sites |±1). To be compatible 
with the string order, spins must be arranged with uOdO 
order which is no more hidden. Thus the strong D-term 
selects the spin states among the states with hidden order 
to realize the explicit Neel order of various periodicity. 

Near the line D ~ SD » J, we can also derive the 
effective Hamiltonian. The spins on the odd-th sites are 
confined into the state |0) and those on the even-th sites 
can take all |±1) and |0) states. The effective coupling 
between the latter spins are obtained within the second 
order perturbation in J and first order in D — SD as, 
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For the model (4), it is known that the large- D phase is 
stable for D cS /J cS > 0.35. For 0.35 > D cS /J cS > -2.0, 
the ground state is on the critical line between the Hal- 
dane and XY phase. For D c s/J c s < —2.0, the ground 
state is on the critical line between the Neel phase and 
the XY phase. 2-5 In the latter two cases, considering 
the continuity to other parts of the phase diagram, we 
may expect that the ground state is actually the Hal- 
dane phase or the Neel phase due to the Ising compo- 
nent of the effective exchange interaction which appear 
as higher order corrections in J. It should be noted that 
this Neel phase is the uOdO phase in the original Hamilto- 
nian (1). The large- D-Haldane phase boundary is given 
by Defi/Jeff = 0.35 which gives, 
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The uOdO-Haldane phase boundary is given by 
Defi/Jefi = — 2-0 which gives, 

= -6.0J 2 (7) 



Dl - SD 2 



These lines are plotted by the broken and dotted lines 
in the ground state phase diagram in Fig. 1 along with 
the numerical results obtained in the next section. We 
thus expect that the Haldane phase shrinks but does not 
vanish even in the limit of large Do ~ 5D. 

On the other hand, for large negative D , the ground 
state is the almost perfect udud state for D < SD and 
the almost perfect uOdO state for D > SD. The ground 
state energy for the former case is given by 



-Budud = N{D - J) 



and in the latter phase 
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E u0dQ = -(D Q - SD) 
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Therefore the transition between these two phases should 
take place around SD ~ — D +2 J. Around this line, how- 
ever, the hidden antiferromagnetic order cannot freeze, 
because the correlation between the 2i-th spin and 2i + 2- 
th spin can be ferromagnetic or antifcrromagnetic de- 
pending on whether the 2i + 1-th spin is |0) or |±1). 
Therefore the hidden order remains hidden. Accordingly 
a narrow Haldane phase should survive in the neighbour- 
hood of the line SD ~ — Dq + 2 J for arbitrarily large SD. 

3. Numerical Analysis 

3.1 Haldane- Large-D transition line 

This phase transition is a Gaussian transition. In or- 
der to determine the phase boundary with high accu- 
racy, we use the twisted boundary method of Kitazawa 
and Nomura. 6 ' 7 The Hamiltonian is numerically diag- 
onalized to calculate the two low lying energy levels 
with the twisted boundary condition (Sf I+1 = —Sf, 
S y N+1 = Sf, S Z N+1 = Sf) for N = 8,12 and 16 by 
the Lanczos method. 

For small D and SD, the ground state is the Haldane 
phase with a valence bond solid (VBS) structure. 8 Under 
the twisted boundary condition, the eigenvalue of the 
spin reversal T (Sf — > —S*,S^ — > — Sf) is equal to 
— 1 in this phase. 6 ' 7 As Do increases with positive value, 
a phase transition takes place from the Haldane to the 
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Fig. 2. Energy gap AE on the line SD = —Do + 2J plotted 
against SD. The vertical axis is multiplied by SD to clarify the 
<5D-dcpcndcncc of AE The data for N = 12 and N = 16 over- 
wrap within the size of the symbols. 



Fig. 1. Ground state phase diagram of the Hamiltonian (1). The 
error bars are within the size of the symbols. The solid lines 
are the guide for eye. The broken and dotted lines represent the 
approximate phase boundary given by Eqs. (6) and (7), respec- 
tively. The dash-dotted line is the line SD = —Do + 2 J. 



large-D phase for which T = 1. Because the quantum 
number T of Haldane and large- 1? phases are different for 
twisted boundary condition, the ground state energy of 
these two phases cross each other at the phase boundary. 
For example, if SD is fixed, the two levels cross at the 
finite size critical point D 0c (N). Extrapolation to the 
thermodynamic limit is carried out assuming the formula 
D 0c (N) = £> 0c (oo) + aN- 2 + WV" 4 . 

3.2 Haldane-udud and Haldane-uOdO transition lines 

From symmetry consideration, these transitions are 
Ising type transitions. We therefore employ the phe- 
nomenological renormalization group (PRG) method to 
determine the transition lines. The Hamiltonian is nu- 
merically diagonalized to calculate the lowest energy gap 
AE(N, D , SD) with total S z = in the periodic bound- 
ary condition using the Lanczos algorithm. Both Neel 
type states have S z = and are two fold degenerate in 
the thermodynamic limit. For finite N, this degeneracy is 
lifted and the energy difference between them gives the 
smallest gap AE(N,D ,SD) which decreases exponen- 
tially with N. On the other hand, in the Haldane phase, 
the energy gap AE(N, Do, SD) remains finite in the ther- 
modynamic limit. Thus the product NAE(N,D Q ,SD) 
increases (decreases) with N in the Haldane (udud or 
uOdO) phase. The intersection of NAE(N,D a ,SD) for 
two successive sizes N — Ni and N 2 defines the fi- 
nite size critical point SD C (N\, N 2 ) for fixed Do- Tak- 
ing into account the Ising universality class, these val- 
ues are extrapolated using the formula 5D c (Ni, N 2 ) = 
SD c (oo) + 2C 1 /(N 1 + N 2 ) 3 . 9 ' 10 The same procedure can 
be also carried out interchanging the roles of SD and Do- 
For the Haldane-udud phase boundary, the system sizes 
N = 8, 10, 12, 14 and 16 are used. For the Haldane-uOdO 
phase boundary, only the sizes of multiples of 4 are com- 
patible with the uOdO order. We therefore use the data 



for N = 4, 8, 12 and 16 to determine this phase boundary. 

The numerical phase diagram also shows that there 
exists an intermediate phase between the uOdO phase and 
the udud phase which is continuously connected with 
the Haldane phase. This is consistent with the analytical 
argument in the preceding section. To further confirm 
the presence of this phase, we calculate the energy gap 
AE to the lowest excited state which has total S z = 1 
along the line SD = —Do + 2 J up to SD — 40 J as shown 
in Fig. 2 for N — 8, 12 and 16. It is clear that the energy 
gap does not depend on the system size significantly and 
behaves as SD -1 for large SD. We have also checked that 
the parity T with respect to the spin inversion is always 
— 1 with twisted boundary condition along the line SD = 
—D + 2 J. This confirms the VBS nature of this phase. 

4. Long Period Modulation 

The mechanism of the selection of uOdO-phase given 
in §2 suggests that the similar mechanism also works for 
a variety of modulation patterns of D-terms producing 
a variety of long range Neel order. To confirm this idea, 
we also carried out the numerical calculation with the 
Hamiltonian with period p modulation of the Z?-tcrms. 
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where Si are spin-1 operators. For p = 2, this reduces 
to the Hamiltonian (1) with Do = 0. The Haldane-Neel- 
type transition takes place at SD C ~ 2.69 for p = 4 (es- 
timated from N = 8 and 16), SD C ~ 2.75 for p = 3 
(estimated from N = 6, 12 and 18) and SD C ~ 3.29 for 
p = 2 (estimated from N = A, 8, 12 and 16). 

For SD > SD C , only the spins S p i are alive. The spin- 
spin correlation functions among these alive spins are 
plotted in Fig. 3 for SD = 4. It is clear that the alive 
spins are aligned antiferromagnetically with period p ir- 
respective of the short range correlation expected in the 
absence of <5£>-term. 



4 J. Phys. Soc. Jpn. 



Letter 



Kazuo Hida and Wei Chen 



<S P S pi > 




Fig. 3. Correlation functions between the alive spins (^SpS*^ for 
p = 1(N = 16, circles), p = 2(N = 16, squares), p = 3(JV = 18, 
diamonds) and p = 4(7V = 16, triangles). <5D = 4. 



5. Summary 

The ground state phase diagram of the spin-1 Heisen- 
berg chains with alternating single-site anisotropy is 
determined by analyzing the numerical diagonalization 
data and the perturbation analysis of various limiting 
cases. In addition to the Haldane, large-D and conven- 
tional udud Neel phases, the Neel phase with uOdO spin 
configuration is found. The emergence of these two types 
of Neel phases is explained as a selection of spin con- 
figuration compatible with the spatial modulation of in- 
terims among various spin configurations with hidden an- 
tifcrromagnctic string order in the Haldane phase. 

We find no direct transition between the udud and 
uOdO phase and the intermediate gapful phase is always 
present. This gapful phase is continuously connected 
with the Haldane phase. Between the large- D phase and 
the uOdO phase, we also find the Haldane phase. This ro- 
bustness of the Haldane phase is also understood by the 
above picture. 

The case of longer period D-modulation is also inves- 
tigated to confirm the above picture. It is natural to ex- 



pect various types of Neel order emerge from the Haldane 
state in the presence of various types of spatial modula- 
tion of D. 

So far, it is impossible to observe the string order itself 
experimentally because the spatial spin structure corre- 
sponding to the string order is fluctuating quantum me- 
chanically. However, if the S = 1 chain material with 
spatially modulated single-site anisotropy D whose mag- 
nitude is larger than the exchange coupling J is realized, 
the observation of the Neel order in such material can be 
regarded as the observation of the frozen string order. 

The computation in this work has been done using 
the facilities of the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo, the Information 
Processing Center of Saitama University. The diagonal- 
ization program is based on the TITPACK ver.2 coded 
by H. Nishimori. This work is supported by a research 
grant from a Grant-in- Aid for Scientific Research from 
the Ministry of Education, Science, Sports and Culture 
of Japan. 
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